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the QCD infrared behavior. 
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I. INTRODUCTION 



One of the main theoretical challenges to understand the B meson phenomenology resides in the 
matrix elements calculation of the non-leptonic decay channels. Several methods were developed 
in order to deal with these decays. The first to treat this problem were Bauer, Stech and Wirbel, 
proposing the factorization assumption (FA) and applying it successfully to several decays Q-Q- 
However, this method fails when there are important non-factorizable contributions. Different 
approaches were developed later in order to make progress in these type of calculations. One 
way to proceed, based on the collinear factorization theorem, is the perturbative QCD (pQCD) 
formalism as applied by Brodsky et al. 4], where it is considered that the non-leptonic decays 
of heavy mesons are dominated by hard gluon exchange and the amplitudes are computed based 



on the hard exclusive hadronic scattering analysis developed by Lepage and Brodsky [a, |{|. A 



powerful method to deal with these decays has been proposed by Beneke et al. @-[TJ, which 
results from the union of the FA and pQCD in the collinear approximation, hereafter indicated by 
QCDF (QCD factorization). These authors argue that the transition form factors, that enter in the 
calculation of B meson decays, cannot be obtained through perturbative methods because they are 
dominated by "soft" interactions, which is a problem that the pQCD advocates claim to be solved 
by the suppression due to Sudakov factors. On the other hand, the non-factorizable contributions 
are dominated by hard gluon exchanges. Therefore, in QCDF the terms that are dependent on 
the transition form factors are parameterized as in the FA formalism, and the non-factorizable 
contributions are perturbatively calculated in the collinear approximation. 

Unfortunately, the formalisms based on the collinear approximations fail in the cases where 
there are non-factorizable contributions related to interactions with the spectator quark and for 
annihilation interactions due to the appearance of endpoint divergences, respectively at twist-3 and 
twist-2, indicating the presence of non-perturbative physics and the breaking of the factorization 
hypothesis [llj. In the QCDF formalism these divergences can be parameterized at the cost of 
some uncertainty in the calculation. These contributions are power suppressed in general, and are 
negligible in front of the factorizablc contributions. This means that QCDF can be successfully 
applied by simply neglecting such contributions. However, annihilation diagrams can generate 
strong phases that are relevant for the CP violation. Besides that, there are decays that occur 
only through annihilation diagrams. These arc quite rare decays, and most of them are beyond 
the present experimental limits. It is expected that they will be detected soon, with the advent of 
LHCb (LHC, CERN), or even earlier by the experiment CDF (Tevatron, Fermilab). These are the 



* Present address: Instituto de Fi'sica, Universidade de Sao Paulo, Caixa Postal 66318, 05389-970, Sao Paulo, SP, 
Brazil. 



2 



decays that we expect to be sensitive to the infrared QCD behavior. 

In this work we compute some of the decay channels of B®(B°) and B®(B®) which occur only 
through the annihilation diagrams: B® — > ir + n~ , Z) ± 7r T and B^ — > K + K~ , DfK T . Working in 
the hard scattering collinear approximation of Brodsky and Lepage, we will use non-perturbative 
solutions for the gluon propagator and for the running coupling constant that were obtained with 
the help of the gluonic Schwinger-Dyson equations (SDE). These non-perturbative solutions lead 
to a gluon propagator and coupling constant that are finite in the infrared, and this fact will 
naturally supplant, at leading twist, the divergences that we discussed in the previous paragraph 
by providing a natural infrared cutoff. 

The phenomenology of B mesons decays is known and thoroughly discussed in Ref. [l3j . However, 
the same cannot be said about the use of SDE solutions in the strong interaction phenomenology, 
and for this reason it is interesting to recall some aspects of this approach. The SDE form an 
infinite tower of coupled equations that in order to be solved must be truncated at some order. 
Most of the solutions obtained so far are gauge dependent and differ among themselves. It is not 
surprising that SDE for the QCD propagator and coupling constant lead to different solutions as 
long as they are solved with different truncations and approximations. There is a class of solutions 
for the gluon propagator that goes to D(k) = as the momentum k 2 — > 0, as well as another 
solution that goes to a value different from at the origin of momenta, whose inverse value is 
usually denominated as a dynamical gluon mass. A discussion about the different solutions can be 
found in the introduction of Ref. [141 ] . where the main references are also quoted. Only very recently 
it was claimed that a gauge invariant truncation of the QCD SDE can be obtained systematically 
[l5| . It is important to stress that the infrared finite behavior of the gluon propagator has also been 
confirmed by lattice simulations [16l . |17| , with strong evidence for a gluon propagator with a value 
at the origin of momenta different from zero Accepting the evidences for an infrared finite 

gluon propagator and coupling constant, now we are faced with the problem of how to introduce 
these quantities in phcnomcnological calculations. 

There are several calculations where the SDE solutions for the gluon propagator and coupling 
constant were used in phenomenological problems, see, for instance, Ref. [l9j]-[28[. In all these 
calculations the fact that the gluon propagator and coupling constant are finite is essential for 
the result, in such a way that no good agreement with the experimental data is obtained without 
appealing to such behavior. This happens, for example, in the case of a non-perturbative Pomeron 
model [H, H3| or in a pion form factor calculation [2l[. Note that an infrared finite coupling 
constant also appears in the context of the so called Analytic Perturbation Theory 29]-[32j, and it 
improves considerably the series convergence of QCD pcrturbative calculations. On the other hand, 
we also expect that any infrared finite gluon propagator leads to a freezing of the infrared coupling 
constant [331 ] . meaning that the use of an infrared finite gluon propagator must be accompanied 
by an infrared finite coupling constant. These finite expressions for the coupling constant and 
gluon propagator should be used in actual calculations in the sense of the Dynamical Perturbation 
Theory (DPT) proposed by Pagels and Stokar many years ago (34| . A finite coupling can also be 
interpreted as an effective charge, as performed in several perturbative QCD calculations [35l [36j . 
Finally, the characteristic mass scale of these solutions is at least a factor 0(2) above the QCD 
scale (Aqcd), and for some of the solutions the free zing of the coupling constant happens at a 
relatively small value [2l|,[37j]- According to Brodsky 38, [Hj], the inclusion of such effects, as the 
freezing of the QCD running coupling at low energy scales, could allow to reliably capture the 
non-perturbative QCD effects at an inclusive level. 

The physics of non-leptonic decays of B mesons may furnish another important test in the quest 
of learning about the non-perturbative behavior of the gluon propagator and the running coupling 
constant. One of these non-perturbative SDE solutions (the one of Ref. [40j]) has already been 
used to compute the decay amplitudes of B mesons [H|]-[43j], although the effect of the associated 
infrared finite coupling constant was not taken into account in these studies. In this work we test 
several solutions and show that they lead to different predictions for the branching ratios of B 
meson decays. Of course, there are many physical quantities that can also modify the decay rates, 
but keeping unchanged all the same quantities (like wave functions and others) while varying only 
the coupling constants and gluon propagators we come to the conclusion that a systematic study 
of these decays will help to select a specific QCD infrared behavior. 

This paper is organized as following: In Sec. 2 we present the basic approach to study the 
non-leptonic B meson decays. In Sec. 3 we show the annihilation amplitudes for the different B 
decays. In Sec. 4 we discuss the different non-perturbative solutions for the gluon propagator and 
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the running coupling constant. Section 5 contains our numerical results and Section 6 is devoted 
to our conclusions. 



II. HAMILTONIAN AND AMPLITUDES FOR B DECAYS 



Weak decays of B mesons are described by an effective Hamiltonian at a renormalization scale 
/i <C Mw- For a B meson decaying into a final state / the effective Hamiltonian is given by (l3j ]: 

where Gf is the Fermi constant, Vckm are Cabibbo-Kobayashi-Maskawa factors, Qi are the oper- 
ators contributing to the decay and Ci(fj) are the Wilson coefficients. 

The operators depends on the flavor structure of the decay. In this work we will discuss two types 
of B meson decays: the charmed decays characterized by transitions with AB = 1, AC = ±1, and 
charmless decays characterized by AB = 1, AC = 0. 

For charmed decays, there are only contributions from current-current operators: 

Qi = (biCj)v-A.(uirj)v-A, Qi = (6»Ci)v-A(%rj)v-A (2) 

For charmless decays, the operators can be divided as: 

1. Current- current operators: 

Ql = (fi«j)v-A(%&t)v-Aj Q2 = (fiUi)v-A.(Ujbj)y-A.. (3) 

2. QCD penguin operators: 

Q3 = fc&Ov-A ^2{gjQj)v-A ; (4) 

9 

Q5 = (nbj)v-A y^fe-g;)v+A ; 

9 

Qe = (fibj)v-A^2(q~jqi)v+A ■ 



3. Electroweak penguin operators: 

3 



Qs = ^(nbj)v-Az2e q ( 



(r 4 6 4 )y_ j4 ^e 9 (g i g i )v+A ; (5) 



2 



3 

9 

9 

where are color indices, r = d or s, e q is the electric charge in |e| units and (<fi<72)v±A — 
gi7 M (l±75)q 2 - The current-current operators describe the W boson exchange at tree level, whereas 
the penguin operators occur at the 1-loop level and describe flavor changing neutral currents. 
The decay amplitude for a B meson into a final state / is obtained from (fT|) 

A(B^f) = (f\H eS \B) 

= ^|E y CKM^(^)(/|Q t |S)(M), (6) 
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where (f\Qi\B) are the hadronic matrix elements between the initial and final states. Due to 
the existence of long distance interactions between hadrons in non-leptonic decays, the matrix 
elements (Qi) must be obtained through factorization schemes or non-perturbative calculations, 
like lattice simulation, QCD sum rules, chiral perturbation theory, etc. The literature contains a 
long discussion about these matrix elements and their contributions to the different factorization 
methods that we quoted before (FA, QCDF and pQCD). 

In the collinear approximation worked out by Brodsky and Lepage the matrix elements of Eq. |6]) 
are obtained through a convolutions of the hard scattering kernel and the distribution amplitudes 
of the mesons involved in the process. In the case of a two-body nonleptonic annihilation decay B 
— > M1M2, where the final state can be two light mesons or a heavy and a light meson, the matrix 
element of an operator Qi is given by: 

(f\Qi\B)= f d[x] Ti([x]) $Mi 0) $m 3 (y) $b0) , (7) 
Jo 

where [x] = x, y, z are the momentum fractions, $m(^) are the light-cone distribution amplitudes 
for the quark-antiquark states of the mesons, which are non-perturbative functions of the momen- 
tum fraction carried by the partons; Ti is the hard scattering kernel that can be perturbatively 
computed as function of the light-cone momenta of collinear partons. 

The Eq. ([7]) is used in the QCD factorization approach introduced by Beneke, Buchalla, 

Neubert and Sachrajda (BBNS) to compute non-factorizable contributions such as hard spectator 
interactions and annihilation. In order to have an actual factorization, the divergences originated 
in the non-factorizable contributions should be absorbed by the distribution amplitudes. However, 
there are endpoint divergences that break down the factorization. These endpoint divergences are 
certain to appear in the kernel calculation of annihilation contributions. The introduction of a 
cutoff is then necessary, and this is performed by BBNS through the following parameterization 

/-=lrx^(l + ^), 0<g<l. (8) 

These are the cutoffs that we referred to in the introduction. This is an inconsistency of the 
method that introduces a large uncertainty in the calculation, in such a way that we obtain only 
an estimate of the amplitudes. 

The end-point singularities that we discussed above also appear in perturbative QCD calculations 
of B meson transition form factors. This type of infrared problem can be handled in the same way 
as will be discussed in the next sections. This problem was already studied by some of us in the 
cases of the 7 — * ttq transition and n meson form factors [U [23| . The only difference in the B 
meson case is that the effect of dressed gluons will not be as pronounced as it is for the ir meson 
case, due to the large mass of the B meson that now appears in the calculation. 



III. AMPLITUDES FOR NON-LEPTONIC ANNIHILATION DECAYS 

In this work we will compute non-leptonic B decays into two pseudoscalar mesons occurring only 
through annihilation diagrams. A typical decay is shown in Fig.[TJ where the bottom quark decays 
via W exchange or penguin processes, and a gluon that is emitted from any one of the quarks 
involved in the process creates a quark-antiquark pair in the final state, which will be a leading 
order process in a s . The channels we will analyze with these characteristics are: 

• Charmless channels : B° s — > 7r + 7r~; B® — > K + K~; 

Charmed channels: B° d -> DfK T ; B° -> D ± tt = f. 



To compute the annihilation amplitudes we use the QCDF approach, where the matrix elements 
are calculated through the Eq. (0. For the distribution amplitudes we adopt the expressions 



found on Ref. at leading twist order. The decay amplitude A(B — > M1M2) is obtained from 



Eq. ([6]) and (O, and the branching ratio is then calculated from 



(9) 



5 



where tb is the B meson lifetime and p c is the momentum of the final state particles with masses 
mi and mi in the B meson rest frame, given by 



_ V( m B ~ ( m i + m 2 ) 2 )(m B - (m! - to 2 ) 2 )) , , 

As mentioned before, these amplitudes contain divergences when dealt with in the collinear 
approximation, as in the QCDF approach. Such divergences can be eliminated in the pQCD 
method when the transverse momenta of the partons are taken into account |44|-[47|. However, 
in this work we will make use of non-perturbative gluon propagators, obtained with the help of 
SDE. These gluon propagators are finite in the infrared, and so they provide a natural cutoff to 
the singularities that appear in the collinear approximation. This kind of calculation has already 
been performed by Yang et al. [Ilj]-[43| for a specific SDE solution, and our intention is to make 
a more general discussion including other solutions, and verify how the results are dependent on 
the particular choices of SDE solutions. We also introduce the effect of the infrared finite coupling 
constant associated to each SDE solution, which was not considered in previous works. 



A. Charmless channels 



There are two charmless decay channels whose final states are light pseudoscalar mesons that 
occur through pure annihilation diagrams: B® — > tt + ir~ and B® — > K + K~ . These channels 
happen through the elementary processes bs — > uu and bd — > uu, respectively, and the creation of 
a dd pair by a gluon, as is shown in the Figs. [2^ and^. These decays receive contributions from 
the W exchange operator Qi and the penguin operators Q4, Qq, Qs and Q10, given in the Eqs. 
©-([S]), with r — d and r — s for the decays of the B^ and B® mesons, respectively. The operators 
Q2, Qa and Qio have a (V — A) ® (V — A) structure while the operators Qe and Qs are of the type 
{V - A)®{V + A). 

We compute the contribution of each operator for the diagrams of Fig. Q] with the help of 
Eq.©. The diagrams that contribute for the amplitude are shown in Figs. [T^ and[T]3, since the 
contributions from the diagrams Figs.JTJ: and[TJi cancel among themselves. The full amplitude for 
these decays is given by 

A(B — ► /) = —7=fBfM 1 fM 2 'X--^- (V ub V* r Arrcc — Vtb V~t*r ^Penguin), 

with r = d and r = s for the decays of the B% and B® mesons, respectively. The tree level and 




FIG. 1: Feynman diagrams contributing to the amplitude of annihilation decay channels. 



G 




FIG. 2: Quark contents of the decay channels: a) B° — > ir + n , b) B% — > K + K 



penguin amplitudes are 



A 



Tree 



C 2 A X 



and 



-4,-.,,^,,,, = [ 2C 4 + ^) Ai+ (2C 6 + ^ ) T>. 



(11) 



(12) 



The functions A\ and A2 are the contributions from operators of the type (V — A) ® (V — A) and 
(V — A) (g> (V + A), respectively, and are given by 



At 



d[X] (j> Ml (x) 4>m 2 (y) 4>b(z) D g (x, y) [a s (fi h )(y - z)D q (x, y, z) + 
+a s (/i)(l -x) D b (x,y,z)], 



(13) 



d[X],<f> Ml (x) (f>M 2 (y)<pB(z) D g (x, y) [a s (fj, h ) (x - z) D q (x, y, z) + 
-a s (fi)yD b (x,y,z)] , 



(14) 



where d[X] = dx dy dz, the quark propagators D b and D q (q = d and s for the and B® mesons, 
respectively) are given by: 

D b X i x ^y^ z ) = (1 ~ ( x - z )( E - y))> 

D^fay,*) = (x~-z)(y-z), 

and D g (x, y) is the gluon propagator, whose the perturbative expression is -D" 1 = k q = y (1 — x). 

The scales adopted in the calculations are \i = mi, for the contribution from the diagram shown 
in Fig. [T^. where the gluon is emitted from the b quark, and fx^ = y/Ajjrib for the contribution 
shown in Fig. |T|b where the gluon is emitted from the "spectator" quark, with = 500 MeV 11 1 . 
When the effective scale \i < m& we use Aqcd = 225 MeV and when \i = nib we use Aqcd = 300 
MeV, and the difference is to match the values of the coupling constant due to different quarks 
thresholds. At the Eq. ((TTJ), the scales of the Wilson coefficients for each term are equal to the 
ones of the coupling constants. 



B. Charmed channels 



The charmed decay channels occur through W exchange processes, and there is no contribution 
from penguin diagrams. The operators that contribute for the effective Hamiltonian are given by 
the current-current operators Qi and Q2 (Eq. @). The Fig. [3^ and[3]D shows the quark diagrams 
for the decay channels — > D^K , which occur through the elementary processes bd — > cu (or 
bu — > cd). The quark diagrams for the decay channel B® — > are shown in the Fig. [3J; and 

[3]i, occurring though the elementary process bs — > uc. 
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c) _ d) 




FIG. 3: Quarks contents of the decay channels: a) -> D~K + , b) B° D+K~, c) B° -> D~tt + , d) 



The full amplitude for the general processes B° — ► Z) T Af ± (with B = Bj, and -B s , D,M = 
D s , K, D, 7r) will be calculated from 

.4(5° - ^M±) = ^Lf B f D f M7T ^V C KMA B a^ D + M± , (15) 

The CKM factors for each decay channel are: Vckm — V* b V u d for the D~K + channel; V* b V c( i 
for D+K~; V us V* b for D~n + and 14 s V^ b for D + ir~. The functions ^bo^dta/* are given by the 
following integrals: 



A b o^d-m+ = J dxdydz4> D (y)4>M(x)D g (x,y)^a s (jj,f)^Cx(fj,f)+ o j 

x [-(l-r 2 )((l-2r 2 -(l-r 2 )x)D c (x)-yD u (y))] + <f> B (z) x 

x [aM C 2 {fi h ) (1 - r 2 ) ((1 - r 2 ) (1 - x) + r 2 (y - z)) y, z) + 

- a s (fi)C 2 (fi)(l-r 2 )((l + r 2 )(y + z)-r 2 )D b (x,y,z)]X, (16) 



and 



/■l { ( C (ut) 
A B° d ^D+M- = J dxdydz (f>D(y) 4>m(x) D g (x, y)<.a s ((j,f) yCi(fif) H 

x [(l-r 2 )((l-2r 2 - (I - r 2 ) x) D c (x) - y D u {y))] +$ B (z)x + 
+ [a s (fi h )C 2 {fi h )(l - r*) [y - z) D d (x, y, z) - a s ((i)C 2 (n)(l - r 2 ) x 

x ((l-r 2 )(l-x)-r 2 +r 2 (y + z))D b (x,y,z)]X. (17) 



The functions D b (x,y, z), Dd(x,y, z), D u (y) and D c = D c (x) are the quarks propagators, given 
by: 

D^\x,y,z) = y + z+{l-y-z){l-x){l-r 2 )- 
D-\x) = {l-x){l-r 2 )- 

D-\y) = y(i-r 2 ). 

The function D g {x,y) is the gluon propagator, whose perturbative expression is D~ x = k 2 = 
y (1 — x) (1 — r 2 ). Note that the integrals in the former equations (and also in Eqs. (fT5)) and 



8 



(|14p) are dimensionless after the simplification of the dependence in both numerator and 
denominator. 

The scales adopted in the calculations are the same as the charmless case: /i = rrif, for the 
contribution of the diagram where the gluon is emitted from the b quark, and /ih = \J A^mft for 
the contribution where the gluon is emitted from the "spectator" quark, with = 500 MeV. For 
the diagrams of Fig. [TJ; andQJl where the gluon is emitted from the quarks in the final state, we 
use the scale (if — fi/2 — nib/2. The QCD scales are also the same as the charmless case. 



IV. NON-PERTURBATIVE GLUON PROPAGATOR AND COUPLING CONSTANT 



To compute the amplitudes discussed in the previous section, instead of the perturbative ex- 
pression D g — (fc 2 ) -1 which is IR divergent, we will use for the gluon propagator D g (x,y) that 
appears in Equations (fT3|) , (fT4|) . (| 16[) and (fT7|) different solutions of the gluonic SDE obtained in 
the literature. Each one of these solutions is also associated to an infrared finite coupling constant 
that enters in the calculation. This coupling constant is evaluated at the scale /x, which is the 
average energy of the gluon emission, where we may still have a small difference between the dif- 
ferent solutions for the non-perturbative couplings. The propagators have the interesting property 
of being infrared finite eliminating the endpoint divergences, and at high energies they match with 
their perturbative counterpart. Note that, as discussed in Ref. [3], in these calculations we are 
integrating over the gluon propagator in a large region of the phase space, and this integration 
is weighted by different distribution functions, therefore the result will depend non-trivially on 
the formal expression of the propagator, as well as it will depend on different values of the non- 
perturbative coupling constant. Even with the large mass scale involved in the problem (the B 
meson mass) we may expect a signal of the infrared behavior of these quantities. 

Notice that the behavior of the coupling constant and the propagators are intimately connected. 
As shown by Cornall [4(| the product g 2 D(q 2 ) is constant, in a more detailed discussion this can 
also be verified in the recent work by Aguilar and Papavassiliou [HI, [52, 65]. Furthermore Alkofer 



et al. in a series of papers [53|-[55j have advocated that a s (q 2 ) = {g 2 / A^Zoiq^Z^q 2 ) (the Z' s 



are the gluon and ghost propagators dressing). Such type of relation is expected in the dynamical 
mass generation mechanism and can be traced back to the Slavnov- Taylor or Ward- type identities. 
Therefore the determination of the actual infrared behavior imply in the associated effect of the 
coupling constant and propagator! 

It could be asked if there is double counting of non-perturbative dynamics between QCDF and the 
employment of infrared finite gluon propagators, since according to the factorization theorem, non- 
perturbative dynamics in B meson decay processes has been absorbed into distribution amplitudes. 
We stress that this is not the case. The effect of the infrared finite gluon propagators affects 
only the hard part of the hadronic matrix elements. The introduction of an infrared finite gluon 
propagator means that as we lower the energy we cannot neglect anymore the gluon dressing 
that is predicted by the SDE. Therefore dressed gluons provide the natural cutoff for the theory, 
and should be considered in the calculation as anticipated many years ago in the proposal of a 
dynamical perturbation theory by Pagels and Stokar [34 ]. Only the hard scattering kernels (like 
T 1 and T 11 in Eq.(6)) will be affected, and the distribution amplitudes should naturally contain 
the typical non-perturbative dynamics of the processes involving meson bound states. 

The usual treatment of B meson decays involves consistent expansions in a s and in 1/rnb- Such 
level of consistency has not been achieved yet in the case of solutions of SDE for the gluon and 
fermion propagators. First, only quite recently it was developed a systematic gauge invariant 
approximation and truncation method for the gluonic Schwinger-Dyson equation (see the work of 
Ref. [15]). This method is consistent with an a s expansion. Secondly, to also have an expansion 
consistent with a 1/wf, expansion for B meson decays it will be necessary to solve the coupled 
Schwinger-Dyson equations for the gluon and fermion propagators, and this is quite away from the 
present status of the theory. The only point that can be guessed at the moment is that the results 
for infrared finite gluon propagator will be only slightly changed by fermionic effects as far as the 
number of quarks remains equal to the known one. 

In this section we present the different SDE solutions for the coupling constant and gluon prop- 
agator (D(q 2 )). They were obtained as the result of different approximations made to solve the 
SDE for pure gauge QCD. There are two classes of solutions according to the behavior of the gluon 
propagator at the origin of the momenta, and both lead to an infrared frozen coupling constant: 
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(a) an infrared finite gluon propagator which is different from zero at the origin of momenta, and 

(b) an infrared finite gluon propagator that goes to zero at the origin of momenta. Although the 
lattice simulations results definitely point out to an infrared finite propagator, there is no agree- 
ment about its behavior at q 2 = 0. A quite recent lattice argument favors a non zero value for 
D(0) M- 

The SDE solutions for the gluon propagator were obtained in different gauges as well as in 
a gauge independent approach. Here we just assume their validity for any gauge and write the 
propagator as 



(18) 



with r) = in the Feynman gauge and r) = 1 in the Landau Gauge. In the following we briefly 
describe some of the solutions found in the literature. 



A. Infrared finite propagators with D(0) ^ 
1. Cornwall solution 40 

This solution was obtained by Cornwall many years ago and predicts a running coupling con- 
stant and gluon propagator which are infrared finite and non-null at the origin. It also predicts 
the existence of a dynamical gluon mass which is responsible for the infrared behavior of the 
aforementioned quantities. This dynamical gluon mass has a dependence on the momentum given 
by 



m(¥) 



(19) 



where m g is the gluon mass scale and A = Aqcd is the QCD scale. The running coupling constant 
is 



la/ 2\ 

a s [Q ) = 



4tt 



A)ln( 



9 2 +4M2 (q 2 ) 

TP 



(20) 



with /?q = 11 — f 71 /! an d n f ls the number of active quark flavors at a given scale. Although the 
SDE were solved for the pure gauge theory it is usually assumed that their solutions will not be 
strongly affected by the introduction of fermions as long as nf is not too large [H, HJ ■ The gluon 
propagator is equal to 



Du(q 2 ) = 



[q2+M$(q2)]f3 g 2 hi 



,2 +4M 2( g 2) 

A 2 



(21) 



where the parameter q 2 i s the strong coupling given by a s = g 2 /4ir. The dynamical gluon mass 



m g has typical values [19(-[22j, |40(: 



m g = 500± 200MeV 



(22) 



2. Fit of Ref. 50 



Another solution with a dynamical gluon mass was found in Ref. (50j |. This solution was fitted 
by the simplest way to have a dynamically massive gluon, whose mass obeys the standard OPE 
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behavior at high energy (m 2 (q 2 ) oc m^/q 2 ) 

Ala(<Z 2 ) 



M 2 {q 2 Y 



with mo being the gluon mass scale, and M 2 (q 2 ) the dynamical gluon mass given by 



M 2 {q 2 ^ 



q 2 + m 2 , ' 



(23) 



(24) 



The coupling constant is similar to the previous solution, with M g (q 2 ) being replaced by M 2 (q 2 ). 



3. Aguilar-Papavassiliou solution 51 ' 5 

The SDE solution obtained by Aguilar and Papavassiliou results from a quite detailed analysis 
of the gluon propagator using the pinch technique, which maintains the desirable property of 
transversality of the propagator. This solution is represented by 



Aic(<? 2 



l 



q 2 + m 2 (q 2 ) 



The dynamical mass is: 

2/ 2 
to (q 



q 2 + ml 



( P m o \ 



72-1 



with 72 > 1. The running coupling constant is given by 



A) In 



q 2 + f(q 2 ,m 2 (q 2 )) 
A 2 



where the function /(g , to (q )) is given by a power law expression 



f(q 2 ,m 2 (q 2 )) = Pl m 2 (q 2 ) + p 2 - 



m 4 (q 2 ) 

1 2/„2 



q 2 + m 2 (q 2 ) ' ^ 3 [q 2 + m 2 (q 2 )] 



m 6 (q 2 ) 



(25) 



(26) 



(27) 



(28) 



The parameters p, pi, P2, P3 and 72 are fitted numerically. For A = 300 MeV and too = 250 MeV the 
numerical values of the parameters are the following: p = 2.47, 72 = 1.76 p\ — 5.615, p 2 = —8.523, 
p 3 = 3.584; and for to = 600 MeV the parameters are: p = 1.234, 7 = 1.64, p x = 2.894, 
p 2 = -4.534, p 3 = 2.043. 

The Fig. 2] shown the behavior of the propagators and coupling constants of the three solutions 
presented in this subsection. 



B. Infrared finite propagators with D(0) = 
1. Alkofer et al. solution 53 

This solution, obtained in the Landau gauge, has a different qualitative behavior relative to 
the ones we have presented so far. Although this solution also predicts a freezing of the coupling 
constant, and the solution for the gluon propagator is also infrared finite, the propagator vanishes 
at the origin of momenta. The behavior of the running coupling constant is given by 



liar 2\ 
"s (<? ) 



(*a(0) 



and the propagator is equal to 



ln[e + ai(q 2 ) a * +&i(<Z 2 ) fc2 ] 
Z(q 2 ) 



q 



(29) 



(30) 
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FIG. 4: Coupling constants and propagators for the three solutions with infrared finite propagators and 
D(0) / 0, with in g — m = 500 MeV: (la) Cornwall solution, (lb) solution of Ref. 50, (Ic) Aguilar- 
Papavassiliou solution. 



where Z(q 2 ) and R(q 2 ) are fitted as 



and 



R{q) - i+c( q *r+d( q ^- (32) 



The constants are parameters obtained in the fitting of the SDE numerical solution: 
a A (0) = 2.972, a sA (^i) = 0.9676, a x = 5.292 GcV~ 2a2 , 
h = 0.034 GeV~ 2b2 , a 2 = 2.324, b 2 = 3.169, 
k = 0.5953, <5 = -9/4, c = 1.8934 GeV~ 2K , 
d = 4.6944 GcV- 4k . 
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FIG. 5: Coupling constants and propagators for the two solutions with infrared null propagators: (Ha) 
Alkofer et al. solution, (lib) solution of Ref. 51. 



2. Fit of Re}. 54 



A more recent fit obtained by Alkofer et al. [5J, |55| , resulting from different approximations for 
the SDE, leads to the following expression for the function Z(q 2 ) that appears in the equation (f5D| 

^(g^J^-) (a Rt (q 2 ))-\ (33) 



where 7 = (— 13Nc + 4nf) / '(22Nc — 4ny). The the parameter w comes from the normalization 
of the function Z{p 2 ) as Z(p 2 = lGeV) = 1, giving w = 1.32 for A = 300 MeV and w = 1.4 for 
A = 225 MeV. The propagator for this solution is then 

Du b ( q 2 ) = , (34) 

q 2 

and the running coupling constant is given by 

a S (o) , 4tt q 2 ( 1 1 



"fit 



l + q 2 /A 2 /3 Q q 2 + A 2 \lll{q 2 /A 2 ) q 2 / A 



with p = (11N C - 2n/)/3. 

The Fig. \E\ shows the behavior of the propagators and coupling constants of the two solutions 
presented in this subsection. 
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V. NUMERICAL RESULTS 



In the numerical calculations we use the asymptotic expression for the wave functions of the 
light mesons, pions and kaons, 

**,k(x) = 6x(l-x), (35) 

while for the D mesons we use 

$ D (x) = 6x(l-x)[l + a D (l-2x)], (36) 

with do — 0.8 for Z? ± and a,£> — 0.3 for the Df meson [H, We neglected the momentum 
of the "spectator" quark in the calculation, i.e. we assume x, y S> z, and in this approximation 
the distribution amplitude for the B meson is given by a delta function. Some results will also be 
presented in the case where the wave functions of light mesons are represented by an expansion in 
Gegenbauer polynomials. The Wilson coefficients are computed using the equations given in the 
appendices of Ref. [H| . We do not include strong phases in the calculation of the integrals, i.e. our 
amplitudes will be real apart from CKM phases. It is known that the strong phases generated by the 
annihilation amplitudes are crucial for predicting CP asymmetries. For example, in the Ref. [43| 
the strong phases of annihilation amplitudes are generated by applying the Cutkosky rules for 
the quark propagators. However, the asymmetries are obtained from the ratio of amplitudes and 
the CP asymmetry parameters are not sensitive to the cutoff effect (or the infrared finite gluon 
propagator). Moreover the main theoretical errors originate from the CKM parameters (43|. 
We also use the following parameters [58J ]: 



• Masses: 



m Bd = 5.28 GeV, m Bs = 5.37 GeV, m D = 1.87 GeV, m Ds = 1.97 GeV, m b = 4.7 
GeV; 



• Decay constants: 



f Bd = 0.200 GeV, f Bt = 0.236 GeV, f D = 0.226 GeV, f D , = 0.241 GeV, f„ 
0.132 GeV, f K = 0.160 GeV; 



• Lifetimes: 



r Bd = 1.54 ps, t Bb = 1.466 ps; 

• CKM parameters: 

A = 0.818, A = 0.2272, p = 0.221, fj = 0.340. 

The TableUshows the branching ratios obtained for each decay channel for the different propaga- 
tors and coupling constants discussed in the previous section. As we shall detail in the conclusions, 
it is possible to see that the non-leptonic B decays resulting from annihilation channels are sensitive 
to the infrared QCD behavior. 



VI. CONCLUSIONS 



We have studied some decay channels of B® (B°) and B^ (-B°) mesons which occur through the 
annihilation diagrams: £>? — » ir + ir~ ~, D ± nT and 5° -> K+K~, DfK T . We have argued that 
infrared finite gluon propagators and running coupling constants obtained as solutions of the QCD 
Schwinger-Dyson equations, may serve as a natural cutoff for the end-point divergences that appear 
in the calculation of these decays. 

We computed several branching ratios for some of the different solutions of gluon propagators 
and coupling constants found in the literature. These different solutions appear due to different 
approximations performed to solve the SDE. Our results are shown in Table fl] We argue that 
a systematic study of B decays may help to give information about the infrared QCD behavior. 
Note that the results of Table U are not trivial, in the sense that the differences in the branching 
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TABLE I: Branching ratios B for decays of the B meson obtained with the different infrared finite gfuon 
propagators and coupling constants, and with m g — mo = 500 MeV and /i = m&: Cornwall (la), Ref. 50 
(lb), Aguilar-Papavassiliou (Ic), Alkofer et al. (11a), Ref. 54 (lib). We use Aqcd = 225 Mev, for /i < mt 
and Aqcd = 300 Mev, for fi > mj. The experimental data for the branching ratios of the decay channels 
are also shown, and the corresponding references. 

£>(Decay channel) Di a Dib Di c Dn a Dub Experimental data Ref. 

B{B° a ~> 7r+7r-)xl0 7 1.08 1.58 1.02 4.30 3.73 <17 [6Q]~ 

< 13.6 [61] 

B(B^K + K')xl0 8 4.92 7.18 4.63 19.62 16.90 9±™±1 [57] 

< 37 [58] 

4±15±8 [59] 

S(B?-»D"7r + )xl0 6 1.03 1.54 0.99 4.40 4.02 
B(B^D+ir-)xW 7 1.28 1.88 1.22 5.31 4.76 

B(B^DJK + )xl0 B 1.34 1.98 1.27 5.61 3.72 4.6ti'i±1.0 [62] 

3.2 ±1.0 ±1.0 [63] 

3.1 ±0.8 [58] 

g(Bg->D+A'~)xl0 8 0.67 0.99 0.64 2.66 1.79 < 1.1 x 10 5 [58] 



ratios come out from the integration over different expressions for the gluon propagators multiplied 
by different wave functions. It is possible to see in Table Q] that the different classes of gluon 
propagators, according to its infrared behavior, lead to branching ratios differing by a factor of 
approximately 2 to 4. The existing experimental data of the decay channels that we have discussed 
is shown in Table HI The experimental data, when confronted with the results of Table [H may 
already be used to claim that some approximations to solve the SDE of pure gauge QCD lead 
to poor predictions for some of the branching ratios (see, for instance, the results related to the 
propagator indicated by Du a ). 

In the calculation of B mesons decays through the factorization theorems, the dominant sources 
of theoretical uncertainties are the mesons wave functions and the scale parameter /z, as has been 
pointed out in Ref. [Til ] . We have analyzed these uncertainties in our calculations for the decay 
channel B° s — » D~K + and our results are shown in Table ITTl and as we shall see in the sequence all 
these uncertainties are less important than the differences originated by the two classes of gluon 
propagators. We can see that the branching ratios are not very sensitive to the shape of the 
wave functions, with a difference of about 2% between the results using the asymptotic function 
and the expansion in Gegenbauer polynomials, which is much below the present experimental 
accuracy. The main source of uncertainty is the scale parameter, with a variation of about 40- 
60%. This strong dependence is expected since the annihilation decay processes are of 0(a s (fj,)) 
in the perturbative expansion. This dependence can be reduced considering higher order terms 
of the perturbative expansion. Once the uncertainties in these quantities become well known we 
certainly will be able to extract information on the gluon mass parameter. It is important to 
notice that the amplitude that we calculate for each specific decay is a convolution of the gluon 
propagator and coupling constant with the wave functions of the final state mesons, therefore the 
result is peaked at different momenta, depending on the meson masses, and help to discriminate 
the form of the gluon propagator and infrared coupling constant. Actually, this is the reason for the 
different sensitivities on the IR cutoff in the decay channels discussed in this work, which will help 
to eliminate the possible uncertainties. It should also be remembered that the next order twist can 
modify the results, but if we have a large confidence level in the values of these many parameters 
or functions, we certainly can test the infrared QCD behavior in these decays. In Table ITI1 it is also 
evident the differences between using SDE solutions of the type I and II for the propagators and 
coupling constants. 

It is clearly necessary to collect much more data to settle the question about the infrared QCD 
behavior, at least in what concerns the infrared behavior of the gluon propagator and coupling 
constant. We know that the comparison with the experimental data will also depend on quantities 
like the distribution functions, but with the LHCb experiment the uncertainty in the experimental 
data will be narrowed, and it will be possible to explore different models until a reasonable physical 
picture of the B meson decays is obtained. 

From the tables shown above we see that the results depend sensitively on the different solutions 
of infrared finite gluon propagators. It must be said that a lot of progress has been done in the 
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TABLE II: Branching ratios for the decay channel B" -> D~ K+ , B(B° d -> £>7iv + ) x 10 5 . We compare 
two different scales, ii — mi, and ii — m,(,/2, and use the asymptotic wave function of the K meson and 
also the wave function expanded in Gegenbauer polynomials (at the scale fi = mj,, with «i(l GeV) = 0.17 
and 02(1 GeV) = 0.2, Ref. 56). We also show the results for different gluon masses. 





Asymptotic 




Gegenbauer 


Gluon propagators 


fi = m b fi 


= m b /2 


fj, = m b 


Di a , m g = 400 MeV 


1.95 


3.25 


2.00 


Dia, m g = 500 MeV 


1.34 


2.09 


1.36 


D Ib , m = 400 MeV 


2.70 


4.89 


2.75 


Di b , m = 500 MeV 


1.98 


3.48 


2.02 


D Ic , m = 500 MeV 


1.84 


3.10 


1.87 


D Ic , m = 600 MeV 


1.40 


2.40 


1.41 


D IIa 


5.61 


16.78 


5.75 


D IIb 


3.72 


7.32 


3.95 



subject of SDE for the gluon propagators recently, this also happened to attract the attention of 
lattice researchers, and we believe that the lattice simulation of the gluon propagator can only be 
fitted nicely by the "massive gluon solution", see for instance Ref. (65[, although this is not an 
unanimous opinion. If with improved lattice results we come to the conclusion that the massive 
gluon propagator is the one chosen by Nature, we verify that the theoretical uncertainty will be 
only due to the poor knowledge of the gluon mass value. In this case we can see from our results 
that the variation of the branching ratios is not so large, and will be smaller as long as this gluon 
mass is better determined through all its phenomenological consequences. 

In order to fully understand the non-leptonic B decays that we discussed in this work we see two 
specific directions of study. One of these is the need of finding a gauge invariant truncation in the 
context of the SDE. The importance of constructing such a truncation scheme is the possibility of 
defining a non-perturbative effective charge for QCD, which constitutes a generalization in a non- 
abelian theory of the QED effective charge as discussed in Ref. [l5| . If this step was accomplished 
in all its glory, it would be possible discard one of the solutions of the infrared sector of QCD on 
the basis of gauge-invariance. The other direction of study is that there should be a systematic use 
of infrared finite gluon propagators and coupling constant in non-leptonic B meson decays. Any B 
meson decay involving a gluon exchange will be affected by the IR cutoff that we discussed here, and 
we intend to study other decays as well as to implement the higher twist effects in this calculation. 
The infrared finite gluon propagator provides a natural cutoff for these phenomenological processes. 
As we have seen in Table U the branching ratios of heavy mesons depend on the infrared QCD 
details, and this will certainly happen when we have light mesons in the final states. The most 
remarkable result of Table U is that with the concept of a dynamically massive gluon it is possible 
to predict the branching ratios of some B decays compatible with the experimental data, without 
the help of any ad hoc cutoff and with the same mass scale that fits the experimental data of many 
other processes. 
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